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ABSTRACT 
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^ ^ Max-Planck-Institut fur Astronomie, Konigstuhl 17, Heidelberg, D-69117, Germany 
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r"! I We examine two positions, ONI and 0N2, within the Ophiuchus cloud LDN 1688 using 

O f observations made with the ISOPHOT instrument aboard the ISO satelHte. The data include 

mid-infrared spectra (~6-l 2 fim) and several photometric bands up to 200 fim. The data probe 
the emission from molecular PAH-type species, transiently-heated Very Small Grains (VSGs), 
and large classical dust grains. 

^ _ We compare the observations to earlier studies, especially those carried out towards an 

isolated translucent cloud in Chamaeleon (Paper I). The spectra towards the two LDN 1688 
positions are very similar to each other, in spite of position ONI having a larger column den- 
sity and probably being subjected to a stronger radiation field. The ratios of the mid-infrared 
features are similar to those found in other diffuse and translucent clouds. Compared to paper 
I, the 7.7/1 1.3 /im band ratios are lower, ~2.0, at both LDN 1688 positions. A continuum 
is detected in the ~10/im region. This is stronger towards the position ONI but still lower 
^|-' ■ than on any of the sightlines in Paper L The far-infrared opacities are higher than for diffuse 

medium. The value of the position 0N2, T2oo/N{H) = 3.9 x IQ-^^cm^/i^, is twice the 
value found for ONI. 

The radiation field of LDN 1688 is dominated by the two embedded B type double stars, 
p Oph AB and HD 147889, with an additional contribution from the Upper Sco OB associa- 
tion. The strong heating is reflected in the high colour temperature, ~24 K, of the large grain 
emission. Radiative transfer modelling confirms a high level of the radiation field and points 
to an increased abundance of PAH grains. However, when the hardening of the radiation field 
caused by the local B-stars is taken into account, the observations can be fitted with almost 
I no change to the standard dust models. However, all the examined models underestimate the 

level of the mid-infrared continuum. 

Key words: Infrared: ISM - ISM: clouds - ISM: molecules - dust, extinction 
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1 INTRODUCTION 



* Based on observations made with ISO, an ESA project with instruments 
funded by ESA member states (especially the PI countries: France, Ger- 
many, the Netherlands and the United Kingdom) and with the participation 
of ISAS and NASA. Herschel is an ESA space observatory with science 
instruments provided by European-led Principal Investigator consortia and 
with important participation from NASA, 
t E-mail: mrawling@nrao.edu 
t Closed down on 3 1 December 2009 



1.1 Background 

The mid- to far-infrared (MIR, FIR) emission spectrum of the 
diffuse interstellar medium (ISM) and interstellar clouds has fre- 
quently been described via a three-component IR dust model (e.g. 
IPuget & Legei]|l989h . Such models typically feature a mixture of 
(i) large aromatic organic ions or molecules producing the so-called 
Unidentified Infrared Bands (UIBs o r UIR bands, also termed In- 
frared Emission Features, or lEFs) (Bakes et al. 2001 and refer- 
ences therein), (ii) a population of transiently-heated Very Small 
Grains, or VSGs, (ISellgrenlll984l) generating mid- to far-IR emis- 
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sion, and (iii) larger 100 — 2000 A) "classical" dust grains in 
thermal equilibrium emitting in the far-IR (A ^ 80/im; e.g. MatMsl 
ll996L lLr& Greenberg 1997). 

InlRawlings et al. (2005) (hereafter Paper 1), we presented an 
analysis of the relative contributions of these three components for 
G 300.2 -16.8, a local, isolated, high-latitude translucent cloud. 
Through a comparison of the relative contributions of the three 
components, it was demonstrated that the energy requirements of 
both the observed emission and optical scattered light could be ac- 
counted for by the incident local interstellar radiation field (ISRF) 
alone. 

At the same time that G 300.2 -16.8 was observed, compara- 
ble observations were taken of other Galactic IR em ission regions 
using the ISOPHOT instrument ("Lemke et al."l996") aboard the In- 
frared Space Observatory (ISO) (Kessler et al. 1996), to allow di- 
rect comparisons to be made between environments with strongly 
differing ambient ISRF strengths. One such region was LDN 1688, 
the main molecular cloud of the p Ophiuchi cloud complex. There 
have been relatively few studies of high radiation field environ- 
ments such as LDN 1688 in which a single dataset has been used 
to compare the mid- to far- IR emission spectra. Besides its wide 
wavelength coverage 3 - 240 /xm, ISOPHOT also had a large num- 
ber of filter bands covering the mid- to far-IR wavelength range 
more uniformly than its successors (e.g. Spitzer, Werner et al .120041 
and Herschel, Pilbratt et al..2010.) . important for specifying the dif- 
ferent dust components. 

1.2 LDN 1688 and its dust emission 

The p Oph star-forming cloud complex has been a target of 
ma ny studies (for a review of th e properties of the region, see 
e.g.lWilking. Gagn e & Alle"nll2008 V and contains numerous newly 
formed stars and embedded YSOs (Bontemps et al. 1996). It is now 
known to contain a star-forrning cl uster with ^100 stars with ages 
< 1 Myr old ( Pa dgett et al.ll2008h . Due to a combination of the 
presence of associated stars and large column densities of gas and 
dust, LDN 1688 exhibits strong multiple-component IR emission. 
Detailed analysis and modell ing of IRAS observations of this re - 
gion have been conducted bv lBernard. Boulanger & Puge3 (1 19931) . 

We present here ISOPHOT data between 6 and 240 /xm, show- 
ing observations of two ON-source positions in LDN 1688 and a 
nearby OFF position (see Fig. [T]). In the spirit of the aforemen- 
tioned three-component dust model, imaging at 12, 22 and 160 pm 
are shown in Fig.[T]to reflect the distributions of the PAH, VSG and 
large, classical grains, respectively. Position ONI is close to the 
central brightness maximum of the cloud containing the embedded 
star cluster. This position was previous ly found to exhibit a high 
IRAS /i2yum//iooyum fatio jBernard. Boulanger & Puget 199 3). The 
second sightline, 0N2, is located in a dense filament at the Northern 
edge of the LDN 1688 cloud and is exposed to the external radia- 
tion from the Upper Scorpius OB association. It exhibited a low 
IRAS Ii2^mllim^m ratio, implying that the local dust in this posi- 
tion has perhaps been subjected to different environmental effects, 
such as additional photoprocessing. 

We present an analysis of the three dust components along 
the two sightlines. Section [2] details the observations and data re- 
duction. Section [3] summarizes the main observing results and de- 
scribes a semi-empirical modelling of the emission across the broad 
ISOPHOT wavelength range. FIR opacities and gas column densi- 
ties are also estimated. In Sect. IH we present a more detailed phys- 
ical modelling of the mid- and far-IR emission based on radiative 
transfer calculations. The discussion of the results in Sect.[5]con- 
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Figure 2. Flux-calibrated, sky- subtracted, averaged PHT-SL spectra for the 
ONI (black line) and ON2 (red line) positions. 
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Figure 3. A comparison of the ISOCAM-CVF spectrum with the 
ISOPHOT spectroscopy and photometry of the ONI position (filled circles). 
The fine solid line is the PHT-S spectrum (with a correction applied to ac- 
count for differences in pixel binning), the dashed line is the ISOCAM-CVF 
spectrum. The two IRAS fluxes are also included for comparison (open cir- 
cles). ISOPHOT data at A ^ 20/im have been colour corrected. 



tains two main issues, (i) the radiation field in LDN 1688 and its 
role for the IR emission, and (ii) the properties of the dust compo- 
nents, especially as compared to the diffuse/translucent sightlines 
in G 300.2 -16.8 (Paper 1). On the basis of several basic assump- 
tions about geometry and extinction. Appendix 15. 1 [ presents a plau- 
sible description of the local ISRF components. 



2 OBSERVATIONS AND DATA REDUCTION 
2.1 Spectroscopic observations 

2x2 raster spectr a were obtained us ing PHT-S spectrophotometry 
mode (AOT 40) (iKlaas et aLlll994l) . The TDTs for these obser- 
vations were 63901741 (OFF), 63901742 (0N2) and 63901743 
(ONI). The data reduction was performed using the ISOPHOT 
Interactive Analysis Program (PIA) Version V9.0 (Gabriel 2000). 
The OFF position spectrum only exhibited Zodiacal Light (ZL) 
emission, and no PAH features were seen. The data were reduced 
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Figure 1. Infrared emission and visual extinction in LDN 1688. The upper two images are WISE jWright. E. L. et alfcoio h maps for the 12-/im and 22-/im 
filter bands. The lower left is the corresponding Herschel PACS ( Poglitsch et al. 2010) Red map, centred on 160/im. These three maps reflect the distribution 
of the PAHs, VSGs and big "classical" grains, respectively. The lower right image is a A( J) map obtained using the 2MASS JHKs colour excesses of stars 
visible through the LDN 1688 central region via the multi-band colour excess method. The pixel size for the A{J) map is 1.4' with a Gaussian FWHM = 4' 
used as a smoothing function for the individual extinction values. All frames shown are 2.5° squared, with a, (5( J2000.0) co-ordinates shown. The positions 
observed by us using the ISOPHOT instrument aboard the ISO spacecraft are indicated by the symbol " x ", where "OFF" is the off-source reference position 
and ONI and ON2 are the on-source positions. The bright local stars providing a significant fraction of the local ISRF are denoted by "+" symbols. The 
contours for the 12-/im filter are as follows: 1000 — 1600 data numbers (DN) in steps of 100. The contours for the 22-/im filter are as follows: 290 — 340 DN 
in steps of 10. The contours for the 160-/im filter are as follows: 0.2 — 1.2 Jy / pixel in steps of 0.2. The contours for the A{J) map are as follows: 0.4 — 5.9 
magnitudes in steps of 0.5. Note the close spatial correlation between the large grains and the extinction map. 



to the AAP level, and then had the OFF position results subtracted. 
Single PHT-SL spectra were obtained by using PIA to perform a 
statistically- weighted average over the four positions. The final 
errors were generated automatically by PIA via the propagation 
of the initial observational errors. It was found that significant 
deviation from the averages of both the continuum and peak levels 
was seen at the OFF (1,1) raster position. This was attributed to 
detector memory effects. The data at this single raster position 
were therefore rejected, and the final averaged spectra were instead 



generated using the three remaining raster positions. The resultant 
spectra are shown in Fig.[2l 

To obtain an independent check on the overall properties 
and calibration of the ISOPHOT dataset, a comparison with 
ISOCAM-CVF spectroscopy data obtained covering the ONI po- 
sition was cond ucted. The red uced and recalibrated ISOCAM dat- 
acube ( Boulan ger et "^1 120051 : Boulanger, private communication) 
was analysed using the Starlink Gaia software. 
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Table 1. Observational ISOPHOT photometry data and accompanying error estimates. 
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Column descriptions: (1) Filter name; (2) Central wavelength, Ac; (3) Filter width, AA; (4) Aperture size and shape; "diam." denotes the diameter of a circular aperture, '135"x 135"' and '181"x 181"' denote the 
size of a square detector pixel; (5) ISO TDTs of observations, on-source and off-source; (6) (ON — OFF) surface brightness, ASu', (7) Mean PIA-propagated internal error, (Jint. as described in Sect. 2.4 of Paper 
1; (8) Mean external error (adopted from Paper 1, Sect. 2.4), crgxT^ as a fraction of ASj^; (9) Absolute accuracy (see Klaas et al. 1994); (10) Filter-to-filter relative calibration accuracy (see Paper 1); (11) Adopted 
combined accuracy. For the ZL-calibrated PHT-P data, this is taken to be ctext + Filter-to-filter relative calibration accuracy. For the PHT-C data, in the absence of an independent (ZL) calibration curve, this is 
taken to be ctext + absolute accuracy. 

^ Values listed for the PI and P2 bands are after ZL calibration. 

^ The filter-to-filter relative calibration accuracy values listed for the C2 data were used to produce the error estimates in parentheses. These were used when estimating uncertainties on FIR temperatures (see Sect. 5.1 
of Paper 1). These smaller values are justified for this purpose, as the two sets of measurements all use the same detector configuration and only differ in the filter used, and are thus not susceptible to cross-detector 
calibration uncertainties. The quoted values are based on an estimate of upper limits on remaining uncertainties due to responsivity variations derived in Paper 1. 

^ Errors are listed in parentheses. 
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2.2 Photometric observations 

The photometric observations were performed in the filters that are 
Usted in Table[T] along with their corresponding PHT-P aperture di- 
aphragm sizes or the PHT-C camera field sizes and the observation 
TDTs. The region corresponding to the ONI position was sampled 
over a square area matching the ISOPHOT aperture in both size and 
position, enabling an averaged spectrum to be extracted. The results 
are shown in Fig. [3] If an ISO-typical error of 17 per cent on the 
absolute flux levels is assumed for both spectroscopy datasets, the 
spectra are found to be in agreement, both in the IR band peaks and 
the continuum regions. The sparse map observing templates (AOTs 
17/18/19: PHT-P and 37/38/39: PHT-C) were used for photometry 
(iKlaas etaLlll"99i) for each filter. 

2.3 Photometry reduction and error analysis 

The p hotometric dat a reduction was performed using PIA Version 
V9.0 (lGabrielll200d) . The data reduction process as described in 
Paper 1 was used: the data were calibrated against the on-board 
calibration source FCSl, and the ZL. This l atter was achi eved via 
the ISOPHOT ZL template spectra of L einert et and the 

colour-corrected monochromatic COBEfDlRBE fluxes at 3.5, 5, 12 
and 25 /xm, under the assumption of blackbody behaviour. Signal 
drift correction and foreground subtraction were also performed as 
in Paper 1. In all of the filter bands, a clear excess signal is seen 
towards the ON positions. As in this earlier work, we list here es- 
timate of internal and external errors. Unlike for the G 300.2 -16.8 
photometry dataset, only single ON and OFF photometry measure- 
ments were obtained, making a statistical average calculation of the 
external errors, ctext, not viable. However, the fractional external 
systematic errors are not expected to be larger than those in Pa- 

i >er I. Since both t hese errors and the systematic calibration errors 
Klaas et al ]|2002") here arise from the same causes, we therefore 
adopt the statistical and systematic error estimates used in Paper 1 
for the corresponding data in this work. 

The adopted combined errors for the average (ON - OFF) sig- 
nal listed in column 11 of Table [T] were obtained by arithmetically 
adding the statistical external errors (column 8) and the filter-to- 
filter systematic (calibration) errors (column 10). 

When estimating the filter-to-filter accuracies at A ^ GO/xm 
an d the absolute acc uracies at all wavelengths, the results of 
the lKlaas etaD(l2002h investigation of ISOPHOT accuracies were 
adopted. The estimated absolute accuracies are given in column 
9 of Table [T] As in Paper 1, although the absolute accuracies for 
A ^ 60/im are ^ 20 — 25 per cent, the filter-to-filter uncertainties 
for a single detector (e.g. C2) may be substantially smalle r due to 
the el imination of most of the sources of error (see del Bur g o et aP 
I200I Sect. 4). Unlike the shorter- wavelength data, the absence of 
any independent calibration curve to aid cross-calibration for the 
A ^ 60/im data dictated that the adopted combined errors for the 
average (ON - OFF) signal listed in column 11 of Table [T] be ob- 
tained by arithmetically adding the statistical external errors (col- 
umn 8) and the absolute errors (column 9). 

2.4 Extinction Mapping of LDN 1688 

Although LDN 1688 is expected to be opaque at V and /, the 
2MASS JHKs data can be used to derive near-infrared (NIR) 
colour excesses of stars visible through LDN 1688, an d hence 
the NIR extinction. As in Lombard!. Lada & AlvesI J2008h and Pa- 
per 1, we have applied the optimized multi-band technique of 



iLombardi & AlvesI bOOlh . adopting A{J)/E{J - H) = 2.66 
an d A(J) / E(H — Ks) — 4.13, following the extinction curve 
of iMathisI ( [T99Qb . The reference area for setting the absolute ex- 
tinction levels was a lO'-radius region at 16^ 0"^ 10.0^ -24° 32' 
0.0"; (J2000.0). The resultant extinction map is shown in the lower- 
right panel of Fig. [T] The extinction for each map pixel was de- 
rived from the individual extinction val ues of stars by apply i ng the 
sigma-clipping smoothing technique of iLombardi & AlvesI (I2OOII) 
and using a Gaussian with FWHM = 4' as a weighting func- 
tion for the individual extinctions. This produced extinctions of 
A{J) = 1.88 ±0.06 at the ONI position and A (J) = 0.98 ±0.06 
at the 0N2 position. The measured extinction at the OFF reference 
position was A{J) = 0.56 ± 0.06. Subtracting this from the ON 
position extinction values yielded A{J) = 1.32 ± 0.08 for (ONI - 
OFF) and A{J) = 0.42 ± 0.08 for (0N2 - OFF). 



3 ISOPHOT RESULTS AND SEMI-EMPIRICAL 
MODELLING 

Figs. [2] -in show the data for the two observed positions as spec- 
tral energy distributions. The MIR spectra exhibit strong 7.7- and 
11.3-/im features. The underlying continuum emission, longward 
of '^lO fim and noted in Paper 1, is also present in LDN 1688. 

As detailed in Paper 1, a semi-empirical fit to the MIR spectra 
is performed using a combination of Cauchy curves and a silicate 
continuum to model the molecular-level IR band emission, together 
with a modified blackbody fit for the thermal emission at longer 
wavelengths by large grains. 

For the IR band emission, the fits are minimization fits to 
the PHT-SL data wherever available. For wavelengths longward of 
the PHT-SL wavelength coverage, the areas under the fitted curve 
were constrained to fit the colour-corrected integrated 12.8 (where 
available) and 16/im ISOPHOT filter fluxes. For both sightlines, 
the resultant fits (shown in Fig.]?]) reproduce the in-band integrated 
fluxes within 20 per cent of the photometrically-measured values. 

The three longest- wavelength photometry points, CI 00, CI 35 
and C200 are fitted using a modified blackbody function of the form 
f^B^v). This function is individually convolved with each of the 
three ISOPHOT filter response curves in turn, and the temperature 
and scaling are adjusted to best fit the photometry. At positions 
ONI and 0N2, the temperatures of ^23.5 K and ^24.8 K, respec- 
tively, were obtained (Table |2]). 

The resulting fits are shown in Fig. |4] and the parameters are 
listed in Table |2] It can be seen from Fig. |4l that in both cases, the 
middle range of the spectrum (^25 — 70 /xm) exhibits additional 
emission that can be explained by the presence of a transiently- 
heated grain component. Physical grain population modelling is 
therefore required, and is described in Sect. 4. We also note that 
the continuum near 10 fim appears to be a much less significant 
contribution than in Paper 1 . 

3.1 Far-Infrared Opacity 

We derive here an estimate of the ratio between the FIR optical 
depth Tem(A) at 200 /xm and the J-band optical extinction A{J), 
T2m/A{J), as well as the average absorption cross section per H- 
nucleon, ^ t{X) /N{H). 

For the case of optically thin emission and an isothermal 
cloud, the observed surface brightness is /(A) = r(A)B(A, Tdust). 
Using the dust temperature values Tdust for ONI and 0N2 as given 
in Tabled we have hence calculated the optical depths. 
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Figure 4. Semi-empirical fitting results. In the left-hand panels, modified blackbody fits of the form v'^Bu(T) to the 100 — 200 /im data points is shown as 
a black line spanning the 40 — 250 /im wavelength range (ONI: upper panels, ON2: lower panels). The PHT-SL spectra are shown in green. Cauchy band 
profiles fitted to the 6.2 — 12.7 /im features are shown in red. A Li & Draine (2001)-style silicate continuum is shown in blue for the 2 — 18 /xm range, and a 
combination of these is shown as a black line for 2 — 18 /xm. The right-hand panels show the details of the emission feature fitting results. 



The total hydrogen column density N(H) is obtained from 
the A{J) values using the the N{H) \s. A{ J) relationship as deter- 
mined by Vuong et al. (2003) from X-ray absorption measurements 
of background stars in LDN 1688, N{H)/A{J) = 5.57 ± 0.35 x 
10^^cm ~^mag~^. The resu lting N(H) and a^oo values are given in 
Table[3l lVuongetalJ (l2003h noted that the ratio for the the general 
Galactic ISM is N{H)/A{J) = 6.4 - 7.8 x 102^cm-2mag-\ 
i.e. some 15 to 40 per cent higher than for LDN 1688 (a 2a effect). 



4 PHYSICAL MODELLING 

As an alternative to the semi-empirical fits presented above, we 
construct here radiative transfer models to explain the observa- 
tions of the positions ONI and 0N2. We use spherically symmet- 
ric models that are intended to describe the average properties of 
these regions. The radiation field consists of the standard ISRF 
(iMathis. Mezger & Panagiairi983h that is scaled, as necessary, to 
represent the higher intensities found in the cloud. We also exam- 
ine the effect of an additional radiation field component that cor- 
responds to the spectrum of a B3V star (Bruzual & Chariot 2003). 
The modelling is completed using the dust model of Li & Draine 
(2001), hereafter referred to as the LD model. The LD model con- 
sists of PAHs, graphite grains, and silicate grains. 

The radiative transfer models were optimised to reproduce the 
observations of dust emission and extinction. In the minimi- 
sation, the photometric measurements were weighted according to 



the error estimates. Below 12 /xm, we use the PHT-SL spectrum. 
Its weight is divided by the number of the frequency points to re- 
duce its overall influence on the fits. We require that the opacity of 
the models (as measured through the centre of the model clouds) 
corresponds to observations. The weight of this constraint was de- 
termined by assuming a 20 per cent relative uncertainty for the ob- 
served A{J) values. 

The LD models were constructed in part similarly to Paper 1. 
The model column density was kept as a free parameter although 
restricted by the A{J) measurements. In the first models, the scal- 
ing of the ISRF, /cisRF, and the relative abundance of the PAHs are 
treated as additional free parameters (upper frames in Figs. [5] and 
[6]). In the best fits, the PAH abundance and the radiation field inten- 
sity are raised far above the default values. For ONI, the best fit is 
obtained with /cisrf = 30 and with the PAH abundance increased 
by a factor of ^5. For the 0N2 position, the corresponding numbers 
are ^27 and ^2. The opacity of the model clouds, as measured by 
A( J), is less than half of the target value for ONI, while for 0N2, 
A{J) is ^10 per cent above the value derived from observations. 

In the second set of models, the ISRF was fixed to the Mathis 
et al. (1983) value but we added another radiation field component 
that corresponds to the spectrum of a star with a spectral class of 
B3V. Fig.|7]compares the SED of the local ISRF with the spectrum 
of a B3V star. The same figure shows the extinction cross-sections 
for the dust components of the LD models. B3V star is clearly an 
efficient source of dust heating, especially in the case of the PAHs. 
In the modelling, the intensity of the B3V component was scaled 
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Table 2. Results of semi-empirical model fitting. Overall errors on the peak 
heights are estimated to be ~ 17 per cent based on a comparison of the 
ONI spectrum with the ISOCAM CVF spectrum, with relative band-to- 
band errors ~ 10 — 15 per cent (as in Paper 1). F(A) values are quoted in 
units of 10~-'^^W cm~'^ i2in~^ST~^ . 



Feature 


Fit parameter 


ONI 


ON2 


6.2-/im 


Central A [/im] 


6.25 


6.28 


band 


Width [/im] 


0.26 


0.24 




F{X) Height 


329.5 


364.4 


7.7-/im 


Central A [/im] 


7.65 


7.62 


band 


Width[/im] 


0.66 


0.53 




F(A) Height 


426.6 


430.0 


8.6-/im 


Central A [/im] 


8.56 


8.59 


band 


Width[/im] 


0.28 


0.31 




F(A) Height 


122.9 


160.0 


11.3-/im 


Central A [/im] 


11.28 


11.20 


band 


Width[/im] 


0.32 


0.22 




F(A) Height 


217.7 


212.5 


12.7-/im 


Central A [/im] 


12.64 




band 


Width[/im] 


0.7 






F(A) Height 


79.1 




7.7/1 1.3-/im band ratio 


2.0 ±0.3 


2.0 ±0.3 


lO-yum continuum level 


29.9 ±2.2 


11.2 ±0.8 



(average of F(A) over 9.5 — 10.5 /im range) 



Classical grains: equilibrium 23.5 24.8 
temperature, T [K] ±f l ±2-2 



Table 3. Visual extinction, optical depth and average absorption cross- 
sections per H-nucleon for ONI and ON2. r2oo is calculated under the 
assumption that rx oq X~'^. 





ONI 


ON2 


A(J) 


1.32 ±0.08 


0.42 ±0.08 


N(H) [1021 cm-2] 


7.4 ±0.7 


2.4 ±0.5 


hooiiy) [MJy sr-i] 


360.0 ± 79 


266.6 ± 59 


T200 [10-4] 


14.8±5-5 




/20oM(J) [MJy sr-i] 


273 ± 62 


635 ± 62 


r20oM(J) [lO-^mag-i] 


11.2±4-2 


22.o±l:° 


^200 = ^200 /N{H) 


2.0 ±0^8 


3.9 ±1.6 


[10-25 cm2 Hnucleon-1] 







with a factor /cbsv, where a value of 1.0 corresponds to a radiation 
field energy equal to that integrated over the full ISRF spectrum 
jMathis. Mezger & Panagia 1983). The free parameters are thus the 
column density of the model and the factor /cbsv (as indicated in 
the lower frames of Figs. [5] and [6]). For ONI, the minimum is 
reached with /cbsv = 18, while for 0N2, the value is /cbsv = 12. For 
the position ONI, the A (J) value is now within 20 per cent of the 
observed value, while for 0N2 the value is ^30 per cent above the 
observed value. The difference in the quality of the fits between the 
first and second set of models is not very significant. Nevertheless, 
the values are lower by ^25 per cent for both ONI and 0N2. 

We also tested the case in which both /cisrf and /cbsv were 
kept as free parameters. The results are very similar to the case in 
which only /cbsv was varied but the value of /cisrf decreases be- 




10.0 100.0 

A (/im) 



Figure 5. Model fits for position ONI using the Li & Draine (2001) dust 
model. The lines show the three emission components and the total inten- 
sity calculated from the model. The red symbols and the red line show the 
photometric observations (without colour correction, i.e., values at the ref- 
erence wavelength assuming a spectrum uljy = constant) and the observed 
PHT-SL spectrum. The blue symbols are the corresponding values that are 
calculated from the model spectrum for the various ISOPHOT filters. The 
free parameters include A^isrf and /cpah in the upper frame and k b3 in 
the lower frame. The total column density is another free parameter but 
restricted by the observed value of A (J) = 1.2. 

low one. However, the /cisrf values are no longer well constrained 
because the radiation field is completely dominated by the B3V 
component. 

5 DISCUSSION 

5.1 The Interstellar Radiation Field in LDN 1688 

Our modelling in Section 4 of the 6 - 200 /xm emission at the 
two sight lines ONI and 0N2 (Figs. \5\ and [6]) indicated that, in 
order to fit the far-IR region, a radiation field enhanced by a fac- 
tor of ^ 20 - 30 is needed if a normal Solar Neighbourhood 
ISRF SED shape (Mathis . Mezger & Panagiat 1983) is assumed. 
The LDN 1688 cloud is, however, known to contain embedded 
early B type stars with an SED strongly increasing towards UV. 

There are two other sources contributing to the ISRF in 
LDN 1688: 

(i) near to mid-IR radiation from the rich embedded clus- 
ter of newly formed stars in th e central molecular and dust core 
dWilking. Gagne & Allenll2008h . and 

(ii) external UV-optical illumination from bright early type 
stars of the Upper Sco OB association, mainly from the northern 
side of the Ophiuchus cloud complex. 
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A (/im) 



Figure 6. Model fits for position ON2 using the Li & Draine (2001) dust 
model (see Fig.[5]for details). The target value for A{J) is 0.39. 




Frequency (Hz) 



Figure 7. SEDs of the local ISRF and B3V stellar emission, together with 
the excitation cross-sections for the Li & Draine (2001) model components. 
The "PAH" line represents the PAHs and the small graphite grains (below 
50 A in size). 

Estimates of the three ISRF components are presented in Ap- 
pendix |a1 

Among the three local components the dominating one 
is the radiation of the two embedded stars p Oph AB and 
HP 147889. The local ISRF also dominates over the general ISRF 
(iMathis. Mezger & Panagiairi983h by one order of magnitude. It is 
difficult to get an estimate for the intervening effective extinction 
in a scattering and clumpy medium between star and the line-of- 
sight dust. Therefore, we are not able to get a quantitative estimate 
better than within a factor of 3 - 4 of the modelled ISRF. Within 



these uncertainties, however, the ISRF model and the observation- 
ally derived values (Sect.© are in qualitative agreement with each 
other. 

5.2 The dust properties 

The ISOPHOT dataset for LDN 1688 describes the IR emission 
spectrum over a large wavelength range, encompassing emission 
from at least three dust grain populations, as also seen in Paper 1 
in the case of the cloud G 300.2 -16.8. The column densities on the 
LDN 1688 lines-of- sight are higher, but by no more than a factor of 
2-3. The radiation field is significantly higher in LDN 1688, how- 
ever, as indicated by the FIR emission that is one order of mag- 
nitude higher than that of the the cloud G 300.2 -16.8 (see also 
iBoulanger et al.|[l996ah . The different environments could be re- 
flected in the properties of the emission from the (possibly ion- 
ized) aromatic molecular species, transiently-heated VSGs and big 
grains. 

5. 2. 1 Mid- infrared emission features. 

The MIR spectra at the positions ONI and 0N2 (see Fig. [3 show 
some differences but no large deviations are observed even com- 
pared to other diffuse and translucent sightlines (the cloud G 300.2 - 
16.8 included). The locations and widths of the PAH features are 
typical of those foun d in the ISM, including regions illuminated 
by nearby stars (see lTielensll2008h . A comparison of the semi- 
empirical fitting results, listed for LDN 1688 in Table [H and the 
diffuse sightlines examined by Kahanpaa etal ] (l2003h also show 
good agreement between the central wavelengths and widths of the 
emission features with the only exception that the width of the 7.7- 
fim feature is smaller in LDN 1688. The 12.7-/im peaks also appear 
to be narrow in LDN 1688 but the feature is outside the ISOPHOT 
spectroscopic range, and is also the most likely feature to be af- 
fected by uncertainties associated with the presence of underlying 
small grain continuum emission. 

In LDN 1688, the 7.7 / 1 1.3 /xm band ratio (ratio of the peaks 
of the features) is ^ 2.0 in both positions. This is lower that the 
values 2.6 - 3.4 found in G 300.2 -16.8. The result can be af- 
fected by details in the fit of the underlying continuum which is 
very strong in the case of G 300.2 -16.8. In G 300.2 -16.8, the con- 
tinuum estimates depend mainly on the 10 fim photometric mea- 
surements whose relative uncertainty was ^^25 per cent. The differ- 
ences between LDN 1688 and G 300.2 -16.8 7.7 / 11.3 /xm ratios 
therefore remain significant. It is well within the variation seen o n 
other Galactic ISM lines-of- sight dChan et al.ll200ll : lTielensll2008h . 
however, and is small compared to the variation between galaxies 
(iGalliano et al.| [2008). The ratios can be attributed to differences in 
the ionization / recombination rates of the carrier species arising 
from the differences in ambient radiation field strengths. 

5.2.2 Mid- infrared continuum. 

As in Paper 1 , our semi-empirical modelling includes a small grain 
component to accommodate the presence of a continuum around 
10 iim. This appears to be correlated with the mid-IR emission at 
16-25 iim of the aromatic emission features. The region around 
10 iim is traced more reliably thanks to the availability of PHT-SL 
spectra and the additional corroboration from the ISOCAM-CVF 
measurement at the position ONI (Boulanger et al. 1996a). The 
continuum level is higher at position ONI than at 0N2. This is 
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most clearly visible in the 10 /im region, but can also be recog- 
nised as higher values around 7.0 /xm and 8.1 /xm between the PAH 
features. In contrast, the total signals (continuum + feature) at the 
peak positions of the PAH features (6.2, 7.7, 8.6, and 11.3 /xm) are 
prac tically the same f or bot h ONI and 0N2. 

iKahanpaa et al.l (l2003h studied a number of Galactic diffuse 
lines-of- sight and, using the average spectrum, derived an upper 
limit of 10 per cent for the 10/im continuum relative to the peak 
flux density of the 7.7 /xm feature. Using the average signal be- 
tween 9.5 and 10.5 /xm as the continuum value for comparison, the 
ratio in LDN 1688 is 9 per cent for ONI and 5 per cent for 0N2. 
There is therefore no indication that the LDN 1688 sightlines sig- 
nificantly differ from the general ISM in this respect. In terms of 
the column density, G 300.2 -16.8 falls between LDN 1688 and the 
case of diffuse sightlines. There, however, the ratios of the 10 /xm 
continuum and the 7.7 /xm peak flux density were ^-^20 per cent or 
above. 

The stronger continuum at the ONI positions of both 
LDN 1688 and (especially) G 300.2 -16.8 must be produced by 
VSGs or possibly even PAH-sized particles. The increased emis- 
sion could result either from a higher abundance of such particles 
or from a greater radiation field intensity. In LDN 1688, the strong 
radiation of early B-type stars combined with the high extinction 
at the ONI position could conceivably change the relative inten- 
sity of the MIR continuum and the UV-excited PAH features. A 
strong effect cannot be expected, however, because both emission 
components cover the same wavelength range. In G 300.2 -16.8, 
the radiation field is normal, and any explanation of the continuum 
must involve changes in the dust properties. 

5.2.3 Comparison of MIR and FIR emission 

The IRAS 12 /xm / 100 /xm and (especially) the 25/xm/lOO/xm 
surface brightness ratios had suggested strong differences between 
the LDN 1688 sightlines, with the 0N 2 position showing lower 
MIR emission. iBoulanger et al.l (Il996ah already estimated that at 
position ONI, ^60 per cent of the IRAS 12 /xm signal could be 
attributed to the continuum. The IRAS data are generally con- 
sistent with the higher spatial resolution ISO data. However, the 
25 /xm IRAS measurement of the 0N2 position is low compared 
to ISOPHOT P2_20 and P2_25 filter measurements and, therefore, 
the difference between ONI and 0N2 is smaller than in IRAS 
data. The difference in the spatial resolution should not play a sig- 
nificant role. The higher resolution WISE and Herschel maps show 
that the 0N2 position coincides with a peak in the surface bright- 
ness (see Fig.[T]) but the values vary by less than 20 per cent within 
the larger area covered by the IRAS beam. 

We can compare the LDN 1688 and G 300.2 -16.8 sight- 
lines quantitatively using the ISOPHOT filters C_100, C_135, and 
C_200 to trace the FIR emission and the filters P_16, P_20, and 
P_25 to trace the MIR emission longward of the main PAH features. 
With the in-band power values (see Table 2 in Paper I), we obtain 
MIR/FIR ratios of 0.39=b0.10, 0.44zb0.09, and 0.15±0.04 for the 
three positions in G 300.2 -16.8. In LDN 1688, the corresponding 
ratios are 0.25 =b 0.06 and 0.19 ± 0.05 for positions ONI and 0N2, 
respectively. Thus, the LDN 1688 sightlines are very similar to the 
position 0N3 in G 300.2 -16.8 regarding the ratio of MIR and FIR 
intensity, in spite of the very different radiation field. In LDN 1688, 
the MIR / FIR ratio is lower in position 0N2 than in position ONI, 
the difference being only at a la level. Using the FIR part of the 
spectra, we derived a large grain temperature of ^24 K in both 
LDN 1688 positions. Based on the SED shape alone, there would 



not seem to be a large difference in the intensity of the heating radi- 
ation field. However, one must also take into account the difference 
in the opacities. For position ONI, the total optical depth was es- 
timated to be twice as high as for the 0N2 sightline. This would 
allow for a larger amount of cooler dust that is shielded from the 
external UV field. 

5.2.4 Dust properties and the radiation field. 

The physical modelling of Sect.|4]helps us to address the question 
of the relative importance of the intensity and spectrum of the heat- 
ing radiation field and the possible variations in dust properties. By 
taking into account dust temperature variations within the model 
cloud, it could also result in slightly different estimates of T200. 
We examined the heating by B stars relative to the illumination 
by a similarly energetic radiation field with an SED correspond- 
ing to the normal ISRF. The spectra of both LDN 1688 sightlines 
could be fitted satisfactorily by assuming a normal ISRF-type SED 
scaled by a factor ^ 30, and by increasing the relative abundance 
of PAH grains. For position ONI, the PAH abundance had to be 
increased almost by a factor of ^ 7, whereas for position 0N2, a 
factor of ^ 2 was sufficient. As discussed above, however, the radi- 
ation field is known to be affected by the high mass stars inside the 
p Oph cloud and by the Upper Scorpius OB association. This jus- 
tified the alternative models in which the radiation field contained 
a significant harder component corresponding to a B3V star. With 
this additional radiation field component, not only could the obser- 
vations be fitted equally well, but in addition, no modification was 
needed in the dust model. This demonstrates the general degener- 
acy between the properties of the radiation field and the dust grain 
abundance. In the case of LDN 1688, it also suggests that the differ- 
ences in the environmental conditions (the radiation field in particu- 
lar) are not necessarily strongly reflected in the relative abundances 
of large grains, VSGs and PAHs. Nevertheless, one can note that 
the models underestimate the intensity between 10-20 /xm, i.e., a 
stronger continuum component is needed. Furthermore, while the 
models with a B3V stellar component and normal PAH abundance 
could fit the shortest wavelengths very well, the observed intensi- 
ties at ONI are underestimated towards the long wavelength end of 
the PHT-SL spectrum. The effect is rather subtle but may be an- 
other indication of the need for a larger VSG contribution. It could 
also suggest that the employed radiation field is already too hard, 
making the continuum too flat below 10 /xm. The actual LDN 1688 
spectra may thus require not only a more intense UV field, but also 
some increase in the abundance of the PAH and VSG grains. In the 
models, this would decrease the /cbsv values and probably improve 
the fits in the FIR regime where the observations appear to be more 
consistent with a lower large grain temperature. 

In Sect.[TTl we derived 200 /xm optical depths of 1.4 x 10~^ 
and 0.9 x 10 ~^ for the ONI and 0N2 positions respectively. The 
radiative transfer models including the B3V star radiation com- 
ponent resulted in similar values, ^ 1.5 x 10~^ for ONI and 
^ 0.8 X 10~^ for 0N2, and thus do not significantly alter the 
previous cr|^o estimates that were (1.9 ± 0.6) x 10~^^ cm^ / H 
for ONI and (4.0 ±1.3) x 10"^^ cm^ / H for 0N2. The opacity 
is thus higher for the lower column density sightline 0N2. Both 
numbers are higher than the values 1.0 x 1Q~^^ (A/2 50 /xm) cm^ / 
H derived for high latitudes teoulanger etal.lll996bh . For the ON2 
position, the value is actually si milar to what has been reported for 
apparently much denser clouds (Hildebrand 1983: 'Beckwith et aP 
1990; Mar tin et al...2012 : Lehtinen et al. 1998; see Table 8 in Pa- 
per I). It can be conjectured (for example) that photoevaporation 
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caused by the strong UV fields has decreased the grain sizes at the 
ONI position. This is not necessarily the case, however, because 
the dust opacities at both positions are still higher than in diffuse 
medium and, because of the lower opacity, 0N2 should be even 
more susceptible to this process. It is more likely that differences 
in volume density play a role. The 0N2 position lies at the centre of 
a filamentary structure and actually coincides with a local column 
density peak that is visible in (for example) the Herschel map 
(see Fig. 1). On the other hand, 0N2 passes through a envelope of 
LDN 1688 that is likely to be more extended along the line of sight. 
In spite of the higher column density, the volume density along the 
ONI sightline could be lower than for 0N2. In this case, the dif- 
ference in opacities could be explained by the growth of grain sizes 
in a dense environment. Both sources have similarly high colour 
temperatures, Td ^ 24 K, and the A{J) values are too low to com- 
pletely shield the dust from the external UV field. Thus, these do 
not seem to be suitable sites for the formation of either ice mantles 
or large dust aggregates. In the case of 0N2, however, the peak ex- 
tinction could be underestimated because A{J) was derived from 
2MASS stars at a relatively low spatial resolution. In the radiative 
transfer models, the dust temperatures were ^20 K or higher, but 
if the density distribution is strongly peaked, the central extinction 
could be several times higher. This would allow the existence of 
a significant amount of dust below 15 K without it having a large 
effect on the SEDs. 



6 CONCLUSIONS 

Our ISOPHOT spectroscopic and spectrophotometric data for two 
positions within LDN 1688 have probed the properties of the main 
dust populations: PAH-like molecular species, transiently-heated 
VSGs and large, classical grains. 

The mid-infrared molecular emission band spectra are not 
markedly different from those seen in other ISM environments, in 
spite of the strong local radiation field. The 7.7 /xm / 11.3 /xm band 
height ratio is ^2.0 at both positions. This is lower than the values 
observed toward G 300.2 -16.8. 

As in Paper 1, a population of molecular species and/or VSGs 
is needed to fit the 10 fim continuum emission. The continuum is 
higher at the position ONI than at the position 0N2. Taking into 
account the high continuum levels in G 300.2 -16.8, this suggests 
that the variation is due to dust properties, rather than a direct effect 
of a high radiation field. 

The strong local radiation field is reflected in the FIR colour 
temperature that is ^24 K for both ONI and 0N2. However, re- 
garding the the MIR/FIR emission ratios, LDN 1688 is still simi- 
lar to the G 300.2 -16.8 0N3 position (Paper 1). The FIR opacity 
(t2ooM(-^) and cr|^o) towards 0N2 is approximately twice as high 
as towards ONI. 

With plausible assumptions of the relative positions of the 
stars and the intervening extinction (see Appendix |A]), we con- 
clude that the stars Rho Oph AB and HD 147889 dominate the 
local ISRF. The nearby Upper Sco OB association has a smaller 
contribution, mainly at the 0N2 position. The total radiation field 
intensity is probably at least one order of magnitude above the value 
in the solar neighbourhood. 

In radiative transfer models, the spectra at ONI and 0N2 can 
be fitted by increasing the local ISRF by a factor of ^20 - 30 and 
the abundance of PAHs by a factor of a few. However, the mod- 
elling also reveals a strong degeneracy between the assumptions of 
the radiation field and the dust properties. The best fits are obtained 



by taking into account the harder radiation originating in the B- 
type stars, almost without any modification of the PAH abundance. 
However, all models underestimate the 10-/im continuum. 

The small differences between the sightlines can be under- 
stood by the higher optical depth of the sightline ONI counter- 
balancing the higher intensity of the radiation field. The higher dust 
opacity towards 0N2, in spite of the lower column density, could 
be related to higher volume density along that line-of- sight. 
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APPENDIX A: THE INTERSTELLAR RADIATION FIELD 
IN LDN 1688 

The local ISRF in LDN 1688 is dominated by two B-type double 
stars. Their distances, spectral types, luminosities, and extinctions 
are: 

p Oph AB: d = lllllo PC. B2V+B2V, 5800 Lq, A{J) = 0.60 
mag, 

HD 147889 AB: d = 118t\l PC, B2IV+B3IV, 7800 Lq, 



A ( J) = 1.45 mag. 

The distances are from the revised HIPPARCOS cata- 
logue of_va^Leeuwe^2po3), spectral types from SIMBAD 
and 



Casassus et al 



"(l2008h . corresponding luminosities from 



de Jager & Nieuwenhuiizen (1987), and extinctions based on 
Carrasco. Strom & Str om (1973); Whittet (1974). 

Using VLB A parallax measurements of two embedded stars , 
DoAr21 and SI, located close to HD 1 47889 j Loinard et all Ew^ 
derived the distance d = I2O.OI4 J P^' which we adopt here as the 
distance to LDN 1688. 

The extinctions through the cloud at the sightlines of the stars 
(cf. Sect. 2.4), A{J) = 0.88 ± 0.05 for p Oph AB and A{J) = 
2.88±0.07 for HD 147889, are clearly larger than the extinctions of 
the stars, p Oph AB and HD 147889 are surrounded by bright and 
exte nded blue reflection nebulae, v dB 105 and vdB 106 (cf. Fig. 
1 of lWilking. Gagne & Allenll2008h . We conclude that both p Oph 
AB and HD 147889 must be embedded in the cloud. Their distances 
are compatible with that of LDN 1688, i.e. 120 pc. Assuming that 
the stars and the dust in ONI and 0N2 are at the same distance of 
120 pc from the Earth (i.e. all are in the plane of sky) we obtain the 
following values for their separations: 

p Oph AB to ONI 1.44 pc, to 0N2 0.97 pc; 

HD147889 to ONI 0.75 pc, to 0N2 3.05 pc. 
The presence and extension of the blue reflection nebulae also 
demonstrate that ONI and 0N2 are exposed to substantial opti- 
cal illumination from p Oph AB. ONI is exposed to illumination 
mainly from HD 147889. 

In order to estimate more quantitatively the contribution to the 
ISRF by the stars, we make use of results bv Habing C1968) for the 
IS RF at A = 100 - 2200 A. For different spectral types, Table 3 
of Habing lists the distance, rioo, from a star at which the 

stellar contribution to the radiation density equals the average So- 
lar neighbourhood ISRF value between A = 1000 - 2200 A (often 
denoted by Go = 1). For p Oph AB and HD 147889, noo = 8.5 
pc and 9.6 pc, respectively. We thus find the following ISRF con- 
tributions in Habing units for the case of no intervening extinction: 

p Oph AB at ONI Go = 35, at 0N2 Go = 77; 

HD147889 at ONI Go = 164, at 0N2 Go = 10. 

For an estimate of the intervening extinction at 2000 A in- 
fluencing the light passing from each of the stars to position ONI 
or 0N2, we adopt, somewhat arbitrarily, half of the line-of- sight 
extinction through the cloud: for p Oph AB A{J) = 0.44 both 
for ONI and 0N2, and for HD 147889 A{J) = 1.44 for 0N2 
but, because of the shorter traversed distance, A{J) = 0.48 (i.e. 
^ 1/3 of the full-distance value) for ONI. These values corre- 
spond to a model in which the cloud's optical depth is trans- 
versely as large as along the line of sight. For the transforma- 
tion from A{J) to A(2000 A) we use the extinction parameter 
value Ry = A(V )/E(B - V) = 4.2, valid for the p Oph 
cloud d Whittet et aD r200D. Part of the stellar Ught lost by extinc- 
tion is returned to the forward direction by scattering. For com- 
plete forward scattering with albedo a, the effective extinction is 
Aef / = A X (1 — a). With a ^ 0.4 and fairly strong forward scat- 
tering with g ^ 0.6 at 2000 A dLi & Drainell200ll) . a substantial 
part of the scattering is not directed forwards. We nevertheless use 
the estimate Aeff = A x (1 — a) = 0.6 x A which is motivated 
by the phenomenon that the penetration of scattered light is sub- 
stantially favoured by a clumpy medium. The effective extinctions 
A(2000A)arethus: 

p Oph AB to ONI and 0N2 1.7 mag; 

HD147889 to ONI 1.8 mag, to 0N2 5.8 mag. 
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With these extinction corrections the the ISRF contributions 
by the two stars obtain the following estimates: 

p Oph AB at ONI Go = 7 , at 0N2 Gq = 16; HD147889 at 
ONI Go = 28, at 0N2 Go = 0.05. 

These va lues refer to the wavelength region 1000 - 2200 A 
jHabing 1968). The corresponding /cssv values referring to the 
total UV-to-near-IR ISRF as defined in Sect. 4, and added up for 
contributions from both stars, are the following: 

kB3v{0m) = 6.2 and kssviOm) = 2.7. 

The uncertainty is at least factor of 2 and is caused mainly 
by the uncertainty in the intervening extinction estimate. With no 
intervening extinction the values would be: 

kB3v(0m) = 26 and A:s3v(ON2) = 11. 

The radiation field estimates obtained in Sect. |4l (Figs. [5] and 
[6l lower panels) are larger by a factor of 3 - 4 when extinction is in- 
cluded. For the (unrealistic) case with no extinction, the agreement 
would be within ^ 30 per cent. 

Besides p Oph AB and HD 147889, other members of the 
Upper Scorpius OB association also contribute to the radiation 
field in LDN 1688. With a distance of 144 pc to the centroid 
(/ ~ 3 50 deg, b ^ 20 deg) of the early type stars in the asso- 
ciation Jde Zeeuw et al.l 1 19991) and using the radiation parameter 
rioo = 67 pc as given in iHabind (1 19681) we find a contribution 
of Go ~ 7 for no intervening extinction. Because the treatment 
with a point source at its centroid is not a good approximation for 
the widely distributed association of stars, we have also calculated 
separately the contributions of the following eight dominating early 
B type stars: C Oph, v Oph , u Sco A, uj Sco, 6 Sco, tt Sco, t Sco, 
A Lib (Suiat haetal1 l2005). The result is Go = 1.9 for ONI and 
Go = 2.2 for 0N2, the largest contribution of Go ^ 0.7 coming 
from the 09 V type star ( Oph. Thus, the contribution by the Up- 
per Scorpius OB association falls clearly below that of p Oph AB 
and HD 147889 but may be relevant for 0N2 which is close to the 
northern surface of LDN 1688 at which most of the radiation from 
the association is i mpinging. 

We note that iLaureiis et al.l (Il995h have estimated the con- 
tribution of Upper Sco OB to the radiation field impinging the 
L134/183 dark cloud complex, located on the northern side (at 
/ ^ 4 deg, b ^ 36 deg, d = 110 pc) of the association but further 
away from it (- 60 pc) than LDN 1688. Using Habing's (1968) 
results they found a contribution at 1000 - 2200 A corresponding 
to Go - 2. 

^The cluster of newly-born stars jWilking. Gagne & AUenl 

l20Q8h embedded in the dense molecular core ^ 10' West of 
HD 147889 is also a potential source of additional ISRF, at least for 
position ONI. Its near- and mid-IR radiation could be sufficient to 
contribute to the heating of the big grains. Using COBE I DIRBE 
maps we determined the mean surface brightness of the cluster at 2 
- 240 p m within a circle of 0.7° radius. An observer located at po- 
sition ONI would see this same surface brightness over about half 
the sky {2tt s teradians). When ni ultiplied with the absorption coef- 
ficient of the lLi & Drains (|2001) dust model, we find that the total 
absorbed energy from this radiation sou rce corresponds to ^ half of 
the energy absorbed from an isotropic iMathis. Mezger & Panagial 
(1983) ISRF (i.e. kisRF = 1.0). We therefore conclude that this 
radiation field component is not important for the dust heating in 
LDN 1688. 



